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ABSTRACT 

Recent optical afterglow observations of gamma-ray bursts indicate a setting and 
distance scale that many relate to star-formation regions. In this paper, we use 
< V/V max > and a set of artificial trigger thresholds to probe several potential GRB 
source density evolutionary scenarios. In particular, we compare a uniform subset 
of BATSE 4B data to cosmological scenarios where GRBs evolve as the comoving 
density, the star formation rate, the QSO rate, and the SN Type Ic rate. Standard 
candle bursts with power-law spectra and a universe without vacuum energy were 
assumed. Our results significantly favor a comoving density model, implying that GRB 
source density evolution is weaker than expected in these evolutionary scenarios. GRB 
density might still follow star-formation rates given proper concurrent GRB luminosity 
evolution, significant beaming, significant error in standard candle assumptions, or 
were a significant modification of star formation rate estimates to occur. 



Subject headings: gamma ray bursts - stars: formation - cosmology: observations 
supernovae: general - quasars: general 



1. Introduction 



Recent afterglow and counterpart observations of gamma-ray bursts (GRBs), particularly the 
redshift-related observations of GRB 970508, GRB 971214 and GRB 980703, have established 
the cosmological origin of gamma-ray bursts (GRBs) (Bloom et al. 1998, Kulkarni et al. 1998, 
Djorgovski et al. 1998, Galama et al. 1998). What sources could create such high-energy outbursts 
with bulk energy of > 10 51 ? The identification of [Oil] emission lines in GRB 970508 and GRB 
980703 strongly suggests that GRBs are associated with vigorous star formation galaxies (Metzger 
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et al. 1997, Djorgovski et al. 1998). The temporal-spatial coincidence of GRB 980425 and 
SN1998bw, a type Ic supernova, has created speculation that some GRBs are related to certain 
types of supernovae (Iwamoto et al. 1998, Woosley et al. 1998), although in the specific case of 
GRB 980425, the gamma-ray flux would then be low. 

Based on these observations, a variety of GRB models have been proposed, including 
'hypernova', the collapse of rapidly rotating massive stars with strong magnetic fields (B ~ 10 15 G) 
in a binary systems (Paczynski, 1997), and the merger of a helium star and black hole (Fryer & 
Woosley 1998). Models that tie GRBs to massive stars implicitly suggest that GRBs are related 
to star formation in some way. 

On the other hand, the acceptable fit of the GRB luminosity distribution to the comoving 
constant rate in Friedman universe has been thought as a consistent result of cosmological origin 
of GRBs prior to the discoveries of after-glows (Mao & Paczynski 1992, Fenimore 1993, Piran 
1992, Wickramasinghe et al. 1993). Previously, Totani (1997) and Marani (1998) have noted that 
a constant comoving rate for GRBs fits measured log N — log P distributions significantly better 
than published rate estimates involving star formation. Examining possible reasons for differences 
between GRB rates and star formation rates (SFRs) estimates is therefore important. 

In this paper, we make use of < V/V max > to trace source density evolution. Although 

< V/V max > was originally created to test homogeneity in the face of multiple detectors and 
detection criteria (Schmidt et al.1988), this statistic remains a sensitive function of homogeneity, 
and hence source density evolution, even when only a single detector is involved. More specifically, 
we construct a set of sub-samples of BATSE GRBs by artificially changing the threshold of BATSE 
in the 4B Catalog. 

Unlike the more traditional method of quantifying evolution by comparing observational and 
theoretical distributions of log N— log P with a x 2 or Kolmogorov-Smirnov test, the < V/V max > 
comparison we employ has no free parameters. The spectral index is fixed by observations, and 
the evolutionary function is well defined by popular cosmologies. The shape of < V/V max > is 
anchored by the value of < V/V max > at the detection limit. By contrast, the shape of theoretical 
log N— log P distributions is not explicitly anchored by the value of P at the detection limit, 
since different combinations of density and peak-flux could give the same limiting value of P for a 
completely different shape of log N— log P. Therefore, previous log N— log P papers have had to 
add extra information to constrain the shape. For example, Horack, Emslie & Hartmann (1995) 
constrained log N— log P by using time-dilation data, Marani (1998) and Mao & Paczynski (1992) 
used < V/V max >, while Krumholz, Thorsett, & Harrison (1998) used a total intrinsic energy. 

Another novel attribute of using < V/V max > to trace density evolution is that plots of 

< y/V max > show relative densities directly. Differences in densities at each redshift will cause 
theoretic < V/V max > values to diverge from observational values. Relative density differences are 
harder to discern from inspection of log N— log P curves. 

To avoid the uncertainties in the SFR, to simplify the theory, and to address more general 
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evolutionary scenarios like SNe and quasar evolution, we adopt broken power-law models for the 
source density evolution function. This approach allows us a uniform framework for addressing 
several different types of source evolution, without losing the character of information inherent in 
estimated source evolution histories. 



2. Method and Results 
2.1. Theory 

For simplicity of calculation, we assume that GRB are standard candles with power law 
photon spectra dL(E) / dE oc E~@ . For a source located at redshift z, the part of the spetrum 
shifted into the detector energy bandwidth (E\ < E < E 2 ) is: 

JEi(i+z) dE 

where Ei and E 2 delimit the band-pass of detection, and L(z) is the intrinsic photon luminosity 
(photons sec" 1 ). Therefore the observed photon flux of a GRB at redshift z is (Weinberg 1972, p. 
421): 

cE 2 (l+z) dL(E) dE 

rW = L W /^( 2 )( 1 + 2 ) = i«|_ii_ (2) 

The additional redshift factor (1+z) is for the correction of time dilation caused by cosmology 
expansion. r(z) is the proper motion distance, in a cosmology without cosmological constant, r(z) 
is 



r(z) = 



H q 2 



zqo + (go -!)(-! + V2g z+1) 
1 + z 



(3) 



Hq is the Hubble constant, qo is the deceleration parameter. The often-used statistical quantity 
< y/V max > can then be written as 

[Z max ( F mln \ 3 / 2 n(z) r 2 (z) d ''i z ) A~ 

V \ _ / f F mtn \ 3 ^\ _ J o V f ) (i+z)^ Vl+H 2 Qlr2{z)/c2 dz ^ 



V m nr \\ F J / f^a. n(z) , r 2 (z) dr{z) r „ 

J JT+Ij^ \j l+ Hl^(z)/c 2 dz dZ 



where £1^ = —k/R^HQ, k = —1,0, 1 for a universe that is respectively open, flat and closed, and 
qo = 0.5(1 — fifc); n(z) is the number of bursts per unit comoving time per unit comoving volume, 
and z max is the maximum redshift at which the dimmest bursts with F = F m i n is detected. 
Therefore 

N(> /} = 4vrn(z) r 2 {z)dr(z) 

Jo l + z ^l + H^ 2 k r 2 {z)/c 2 
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In a Euclidean universe, standard candle sources with a uniform spatial distribution will 
have < V/V m ax >= 0.5. Source evolution can lead to < V/V max > deviating from 0.5 (Schmidt 
et al.1988). In a cosmological model with constant comoving density, the effect of the universe 
expanding is to decrease < V/V m ax > smoothly as maximum redshift z m ax increases. We refer this 
shape as a standard shape (solid line, figure 3). For a source-evolutionary scenario, < V/V max > 
will deviate from the standard shape. Therefore, < V/V max > is an indicator of differences from a 
constant comoving source density. 



The comoving luminosity density of the universe, L2800> which is considered to be a star 
formation indicator, has been estimated from a sample of faint galaxies over the redshift 
< z < 1 in the Canada- France Redshift Survey (Lilly et al 1996). In qo = 0.5, f2 = 1 universe, 
-^2800 ~ (1 + ^) 3 - 9±0 - 75 to z ~ 1 and reaches its peak between 1-2. SFR then decrease with redshift 
up to z ~ 5 (Madau et al 1996). 

A study of giant HII regions in late-type galaxies indicates that the massive stellar progenitors 
of Type Ib/c supernovae are probably in relatively close binary systems (Van Dyk 1996). Type 
Ib/c and hypernovae progenitors have much larger masses than Type la progenitors, and therefore 
a much shorter lifetime. Their counterparts should therefore directly probe the instantaneous star 
formation rate. We will therefore extrapolate the Type Ic supernovae and hypernovae formation 
rates from the Type II event rate and other assumptions (Madau 1998, Sadat et al. 1998). 

Given that GRB source density is proportional to its progenitor formation rate, we adopt a 
broken power-law evolution function as an approximation. 



where uq is the density at z = zq, zq = 1 for SFR, Type Ic supernova and hypernova formation 
rate and zq = 2 for burst sources that follow the formation rate of quasars. 

In §2.4 we will see that the difference between the anticipated and observed < V/V max > 
becomes larger when a\ and | q?2 I increase. In order to avoid overestimating the evolution of GRBs, 
we choose a\ = 3, medium of 3.9 for the SFR (Lilly 1996) and 2 - 2.5 for the neutron star merger 
rate (Totani 1997). This parameterization can be taken as a rough description of star formation, 
Type Ic SN and hypernova rates. a 2 is set to be -1, from the theoretical calculation for Type II 
SN (Madau 1998, Sadat 1998), though the real GRB evolution might be between or stronger than 
these values. For burst sources that follow the formation rate of quasars, we choose a± = 7 and 
a 2 = -4 (Yi 1994). 



2.2. 



Evolutionary Scenarios 
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2.3. Evolution of Observed < V/V max > 

We selected 739 bursts from the BATSE 4B Catalog that triggered on energy channels 2+3. 
Sub-samples were constructed artificially by increasing the trigger threshold C m in in number of 
counts in 64 ms and 1024ms timescale. As C m j„ was decreased incrementally by a given factor, 
bursts with a newly computed C max /C m i n > 1 were selected, C max is the maximum counts in 
64 and 1024 ms timescale, and we choose the larger vaule of C max /C m i n calculated in the two 
timescales. These sub-samples inherited many attributes of BATSE data although with decreasing 
detection rates, and hence the volumes that encompass these GRBs were incrementally increasing. 

< V/V max > for each sub-sample is calculated by < V/V max >=< (C max /kC m i n )~ 3 / 2 >, where 
k > 1. 

2.4. Results 

We set the GRB spectrum index [3 = 1.1 (Mallozzi et al. 1996). The < V/V m a X > evolution 
with redshift in a universe with Q\ = for different evolutionary scenarios are shown in Figures 
1 and 3. < V/V max > for quasar-like object evolution (dash-dot-dot-dot line, at\ = 7, Q2 = —4) 
shows the greatest deviation from the standard shape (solid line, a.\ = «2 = 0, corresponding to 
non-evolving GRBs). Its peak at z = 2 surpasses 0.6 and decreases to at about z ~ 10, while 
the standard shape smoothly declines to near at z ~ 1000. The dashed line {ct\ = 3, = — 1) 
and dash-doted line (qi = 3, ai = 0) in Figure 1 are nearly unchanged at z < 1. The former line 
slope after z > 1 is steeper than standard shape and the later gradually overlaps the standard 
shape. These results show that stronger source evolution of GRBs will enlarge the deviation of 

< V/V max > from the standard shape. The observed < V/V max >= 0.326 for the entire 739 bursts 
in our sample. This occurs at z max ~ 1.76 for GRBs following the comoving standard shape and 
Zma X ~ 3.6 for bursts following quasar-like object evolution. 

In Figure 2, we normalize the theoretical burst number for z max to the number of bursts in 
our sample: 739. The x-axis is the number of bursts detected for different sensitivities, selected 
by the method described in §2.3. The y-axis is the corresponding < V/V max >. The dotted line is 
the observed evolution of < V/V max >. The dash-dotted line represents the 3a error bars obtained 
by the bootstrap method (Efron 1982). The large fluctuation at the left is due to the small 
number in the sample. The solid lines are the predicted distributions for different evolutionary 
scenarios in a flat universe. We can see in Figure 2a that the calculated distribution with comoving 
constant rate can fit the observed data (dotted line) very well, while the evolution of < V/V max > 
tracing Quasar- like object rate (figure 2b) goes well beyond the 3a upper limit. The intermediate 
parameters shown in 2c and 2d are marginally consistent with the 3a error line. 

In Figure 3, we compare the relations of < V/V max > and redshift in a flat (Ofc = 0), open 
(fifc = 0.8) and closed (f2fc = 0) universe for comoving constant rate (dashed lines) and quasar-like 
object rate (dash-dotted lines). For the observational value of < V/V max > (dotted line), they 
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show a slight difference. 

3. Discussion and Conclusions 

The results presented here imply that GRB source evolution is weaker than expected 
for several evolutionary scenarios. Scenarios that imply significantly fewer bright GRBs than 
observed, and significantly more dim GRBs, include density evolutions following SFRs, supernovae, 
hypernovae, and quasars. 

It is difficult to reconcile this result with indications that GRBs do occur in star-forming 
regions. Several possibilities exist: the first is that GRBs, in general, do not occur at redshifts 
significantly in excess of zq, where the event rate reach its peak, and indications that they do 
are misinterpreted. One example of this may be GRB 971214, which has a reported redshift of 
z ~ 3.4 (Kulkarni et al. 1998). To dramatize this, we present Figure 3, a plot of the theoretical 
distribution of < V/V max > for GRBs originating from redshifts greater than 1. The evolutionary 
< V/Vmax > values are clearly much higher than the observed values. 

A second possibility is that GRBs have a very wide intrinsic luminosity function (Horack 
et al. 1995, Krumholz, Thorsett, Harrison 1998), allowing a fraction of them to lie at very high 
redshifts. If this is true, however, the dispersion of the luminosity function must not only be 
greater than the difference between the comoving and evolving rates, but a systematic bias must 
exist to shift the mean. 

A third possibility is that GRBs display a correlated density and luminosity evolution. One 
natural scenario that might allow for both types of evolution is one that admits relativistic bulk 
motion with a high Lorentz factor (i.e. Meszaros & Rees 1997). Resultant beaming creates an 
evolution of < V/V max > with redshift and increases the maximum redshift z max to a larger value 
(Che et al. 1996, Yi 1994). This effect might shield a discrepancy caused by source evolution and 
perhaps adjust the < V/V max > distribution toward measured values (Che et al. 1998). 

Lastly, we cannot ignore potential systematic biases in extrapolated star formation rates 
(Calzetti 1998, Cimatti et al. 1998), driven by incomplete knowledge of dust extinction, reddening, 
and an inability to detect low-luminosity galaxies. Alternatively, GRBs might originate from other 
kinds of sources whose density and/or luminosity evolution is very different from that of known 
sources. 

After much of this work was completed, the cited manuscripts by Lloyd & Petrosian (1998) 
and Krumholz, Thorsett, & Harrison (1998) appeared on the LANL abstract server. 

This research was supported by grants from NASA and the NSF. 
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Fig. 1. — The variation of < V/V max > with the maximum detection redshift z max is shown for 
the 0^ = universe. GRB sources are assumed standard candles with spectral index 1.1. The 
four lines (solid line, dashed line, dash-dotted line and dash-dot-dot-dot line) correspond to four 
evolutionary scenarios: comoving constant rate (a>i = c*2 = 0); Type Ic supernova and hypernova 
scenario (qi = 3, «2 = —1,0), and quasar-like object scenario (a\ = 7, c*2 = —4). The dotted line 
correspond to the observed value < V/V max >= 0.326. 
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Fig. 2. — Comparisons of theoretical distributions of < V/V max > (solid lines) and the observed 
distribution (dotted line) obtained from 4B catalog. Four evolutionary scenarios are depicted: 
comoving constant rate (a); quasar-like object scenario (b); Type Ic supernova and hypernova 
scenarios (c and d) . The dotted line depicts the observed evolution. The dash-dotted lines are for 
the 3a error bars. 
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Fig. 3. — Comparisons of evolution of V/V max with the maximum detection redshift z max in a 
closed (fife = —0.8, left line), flat (fi^ = 0, middle line) and open (£1% = 0.8, right line) universe 
for comoving constant rate ( dashed lines ) and quasar-like object rate (dash-dotted lines ). The 
dotted line correspond to the observed value < V/V ma x >= 0.326. 
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Fig. 4. — A comparison of theoretical distributions for GRBs originating from a redshift greater 
than 1. We set «2 = 0. The distribution for «2 = — 1 is similar. 



